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We calculate the stress tensor, or energy-momentum tensor, form factors of the pion and of axial 
vector mesons in the chiral limit of a hard wall AdS/CFT model of QCD. One (of the two) pion 
gravitational form factors is directly related to the second moment of the pion generalized parton 
distribution, thus providing a sum rule for the latter. As was also the case for vector mesons, both 
the pion and the axial vector mesons appear strikingly more compact measured by the gravitational 
form factor than by the electromagnetic form factor. 
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I. INTRODUCTION 

In this paper we calculate gravitational form factors, 
which are form factors of the stress or energy-momentum 
tensor, in the axial sector using a hard-wall model of 
AdS/QCD. 

The gauge/gravity or AdS/CFT correspondence is 
studied because it offers the possibility of relating non- 
perturbative quantities in theories akin to QCD in 4 di- 
mensions to weakly coupled 5-dimensional gravitational 
theories [l|, |2|. Some applications that particularly in- 
volve mesons are found in @, i, i, S, 0, i, ©, Ha d 
HI GJ, EI [H, HE Ell, El and other works cited therein. 
The mesons studied are mainly vector and scalar mesons 
and topics studied include masses, decay constants, cou- 
pling constants, and electromagnetic form factors. Less 
studied to date are parton distributions, be they ordi- 
nary ones, or transverse momentum dependent ones, or 
generalized parton distributions (GPDs). 

The present authors [19| have studied gravitational 
form factors and the connection to GPDs for vector 
mesons, obtaining sum rules for the GPDs and finding 
that the vector meson radius appeared notably smaller 
when measured from a gravitational form factor than 
from the electromagnetic form factor. We wish to ob- 
tain the corresponding results for the pion, a particle for 
which it may be easier to obtain experimental informa- 
tion about the GPD [2Cj an d for which there is already 
information on the electromagnetic form factor [211 ]. 

Technically, studying the pion is more involved than 
studying the vector mesons because there are additional 
terms in the action involving the chiral fields. The ground 
has been broken by workers who have studied the pion 
electromagnetic form factor (22|, [23|, [2J| . Our work is sim- 
ilar in its basic approach to the latter two references, but 
we have attempted to make the present paper reasonably 
self-contained. Also, as we are studying the axial sector 
to learn about the pion, it requires only a small extra 
effort to also study the axial vector mesons oi, and we 
quote results for these in the body of the paper. We have 
limited ourselves to the chiral limit, where one can obtain 
analytic results for many of the quantities of interest. 

In general, in AdS/CFT there is a correspondence be- 
tween 4-dimcnsional operators 0(x) and fields in the 5- 
dimcnsional bulk <f>(x, z), where z is the fifth coordinate. 



The 4D sources used in the 4D generating function Z\£, 
we will call </>°(x), and 



j\d\ 



cxp I iS, 



~>AD 



d\. 



O{x)^ (x))\. (1) 



The correspondence may be written as 

Z 4D [(f>°}=e lS5D[ ^' ] , 



(2) 



where on the right, S^cz] is the classical action evalu- 
ated for classical solutions <f> c i to the field equations with 
boundary condition 



lim (f> c i{x, z) = z A cj)°(x) . 

z— >0 



(3) 



The constant A depends on the nature of the operator 
O, and is zero in simple cases [25 L 

The original correspondence l| related a strongly- 
coupled, large N c , 4D conformal field theory to a weakly- 
coupled gravity theory on 5D AdS space. In QCD, N c is 
not large, nor is the theory conformal, as evidenced by the 
existence of hadrons with definite mass. Nonetheless, re- 
sults obtained treating N c as large work surprisingly well, 
and one can argue that QCD behaves approximately con- 
formally over wide regions of Q 2 [l(|. We simulate the 
breaking of conformal symmetry, following the so-called 
"bottom- up" approach as implemented in @, Q , by in- 
troducing a sharp cutoff in AdS space at z = zq. The 
unperturbed AdS space metric is 



MjN 



ds = gMNdx dx 



1 



MjN 



;VMNdx dx 



e < z < zq, 
z~ 

(4) 
where tjmn — diag(l, — 1, — 1, — 1, — 1). The z = e wall, 
with e — > understood, corresponds to the UV limit of 
QCD, and the wall located at z = zq = 1/Aqcd sets the 
scale for the breaking of conformal symmetry of QCD 
in the IR region. [Lower case Greek indices will run 
from to 3, and upper case Latin indices will run over 
0,1,2,3,5.] 

The O to operator correspondences of particular in- 
terest here are 0, H, HI 



T^{x) 



A afi (x,z), 



(5) 



where J^ M 



97^75^9 is an axial current with flavor 
index a, A a ^ is an axial field, T M „ is the stress tensor, 
and hfj, v represents variations of the metric tensor, 

g^ v (x,z) = — (ripv + h IJlv (x i z)) (6) 

We will use h^ v in the Randall-Sundrum gauge [2q . 
wherein h^v is transverse and traceless (TT) and also 
satisfies h^ z = h zz = 0. Variations of the metric tensor in 
a TT gauge will only give us the transverse-traceless part 
of the stress tensor. This will determine uniquely one of 
the two gravitational form factors of the pion, the one 
that enters the momentum sum rule, and correspondingly 
4 of the 6 gravitational form factors for spin-1 particles, 
including the two that enter the momentum and angular 
momentum sum rules. 

Relevant details regarding the pion and axial-vector 
mesons, including the wave functions and the two-point 
functions, are worked out in Sec. [Til and Sec. IIIII works 
out the three-point functions, and extracts from them 
the stress tensor matrix elements. Sum rules and stress 
tensor form factor radii are given in Sec. IIVI and some 
conclusions are offered in Sec. [V] 



up to second order is given by [7| 



S A 



I d 5 x^ 



v{zf 



g M " (d M v a - A a M )(d N Tr a - A a N ) 
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LN 



(8) 



where F a KM = d K A a M - d M A a K , with A=(A L - A R )/2. 



B. Equations of Motion 

Using the unperturbed 5-dimensional AdS space met- 
ric, and taking the variation over A a M of equation l[5|). 
one obtains the equations of motion, which are expressed 
in 4D momentum space as 

d z Qfc^x) + jA a vX 9 -^A" uX = 0, (9) 

d z (~ z d z r) +^(* a - a ) = 0, (10) 

-q 2 d z r + ^d z n a = 0, (11) 



II. PION AND AXIAL- VECTOR MESON 
A. AdS/QCD Model 

The action on the 5-dimensional AdS space is Q 



S5D — 



J d 5 Xy/g\ 11 +12 



(7) 



-Tr 



\DX 



1 

471 



'Z\X\ 2 --^{F 2 +F 2 Rl 



This action contains the X field, which corresponds to 
4D operator q~R,qL and, through F^f R 



= 9 M Al R 
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° A L,R 
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m 



iR i4fjJ, also contains the A a L and 
R fields, which correspond to operators J£ = 
qhlpt a q L and J Rfi = qaj^qn respectively. We define 
Am(x, z) = A^ix, z)t a , where the group generators sat- 
isfy Tr(t a t b ) = 8 ab /2. The covariant derivative of the 
X field is given by D M X = d M X - iAjfX + iXAj?. 
Moreover, the X field can be written in exponential form 
as X(x,z) = Xq(z) exp(2zi a 7r a ). Solving the equation 
of motion of Xq(z) 1 one obtains Xq = ^lv(z), where 
v(z) = m q z + (Tz 3 . Using the AdS/CFT prescription and 
the fact that qaqL appears in the mass term of QCD La- 
grangian, parameter m q can be identified as the quark 
mass and parameter a as the quark condensate (qq). 
In this paper, we will discuss only the chiral limit of 
the AdS/QCD model, i.e m^ = case, or equivalently 
m q = 0. 

The axial- vector and pseudoscalar sector of the action 



where the gauge choice A z = has been imposed. The 
<j) field comes from the longitudinal part of A a = A a + 

The fields above can be written conveniently in terms 
of bulk-to-boundary propagators as follows 

A;(q,z) = Af(q) ± A(q,z)+A?(q) n <Kq,z), 

7T a (q,z) = l -^Af(q)n(q,z), 

q z * " 



with Af{q) 



and A?P(q) 



are the Fourier transform of 
the source functions of the 4D axial current operators 
Ja,!_i( x )± an d Ja,/j,( x )i< respectively. 

The n th KK-mode axial-vector meson's wave function, 
denoted by ipA, n (z), is the solution of equation ([9]) with 
q 2 = rn An , and with boundary conditions ip(0) = 
and d z tp(zo) = 0. The normalization of ipA,n is iden- 
tical to that of the p-meson's wave function, and is given 
by J(dz/z)ipA.n{z) 2 = 10- On the other hand, the 
pion's wave function, denoted by <f>(z) or tt(z), is the so- 
lution of the coupled differential equations, i.e., equation 
(fi"0|) and (jll|) . with q 2 = in the limit of massless pion. 
Furthermore, in this limit, the boundary conditions for 
the pion wave functions are </>(0) = 0, 7r(0) = —1 and 
d z 4>{zo) = 0. [A careful analysis is given in ref. [3] for ra^ 
small. There, the UV boundary conditions for the scalar 
wave functions are <fi(0) = and 7r(0) = 0. However, 
the function tt(z) away from z = approaches —1 rather 
quickly. For m^ — * 0, the function n(z) equals —1 in 
essentially the entire slice of 5D AdS space, < z < zq.] 

Using Eq. (fTTj) and the UV boundary conditions, one 
finds 7r(z) = — 1 for all z. Therefore equation (flU)) can be 



rewritten in terms of \&(z) = 4>(z) — ir(z) as 



i Wi-^S* 



(12) 



The solution is given by [24 



*(z) 



zT[2/3] 



I_i(az 3 )-Ii(az 3 )-^ 




where a = gacr/3, with g§ = 2tt as shown in ref 
Note that \&(z) is identical to A(Q,z). Parameter zq = 
1/Aqcd is determined by the experimental value of p- 
meson's mass m p = 775.5 MeV [27J, which corresponds 
to z = 1/(322 MeV) Q. 

For the ai's wave functions, one has to rely on numer- 
ical methods. However, other aspects of the axial-vector 
mesons are analogous to the vector mesons. For instance, 
the bulk-to-boundary propagator of a\ can be written in 
terms of ipA,n as 



A(q,z) 



E 



(jd z <ll)A,n(e))ll>A,n(z) 



(14) 



c.f equation (17) of ref. [H and ref. @,[U|. The bulk- 
to-boundary propagator A(q, z) satisfies A(q, s) = 1 and 
d z A(q,z ) = 1. 



Two-Point Functions 



The completeness relation is given by 
d 3 q 



E 



(27r) 3 2g 



\n(q))(n(q)\ = l 



(15) 



The complete set of states includes \A^(q,X)), the n th 
axial- vector state, as well as |-7r a (<?)), the pion state. 

By applying the completeness relation into 
(0\TJ1(x)±J%(0)_l\Q) and (Q\TJ^{x\ J^(0) n \0), then 
multiplying q 2 — m 2 and taking the limit q 2 — > m 2 , one 
can extract the following quantities from the AdS/QCD 
correspondence [/J 



*U, 



d z lf>A,n(z) 



f = - 

J 7T 



95 z 
d z V(z) 



9\z 



(16) 
(17) 



where Fa ; „ is the decay constant of the n" 1 mode of the 
Oi, and /„• is that of the pion. The latter is obtained in 
the chiral limit q 2 — > m 2 = 0. The decay constants are 
defined by 



(0\J%JO)±\A b n (p,X)) 
<0|Jl,„(0)„|7r b (p)) 



F Ain£fi (j>,\)5 ab , 

U P ,s ab . 



(18) 
(19) 



Equation (|17[) and (| 13[) relate the input parameter <r 
to the pion decay constant 



fl 



3 T(2/3) 
4tt 2 r(l/3) 



. „ 1 i"2(a^) 
2a 2 ^ 3 3 



J-l(azo) 



(20) 



Using the experimental value / ff = 92.4 MeV, we find 
a = 2.28 A 3 QCD , therefore a = (332 MeV) 3 . Conse- 
quently, other observables can be determined, rriA,i = 

1376 MeV and F A f * = 493 MeV. They arc in a 
good agreement with the experimental values rriA.i = 

1230 MeV and F A 'l = 433 MeV. 



III. GRAVITATIONAL FORM FACTORS 

The three-point function that includes the stress tensor 
follows from 



(0\Tjr{x)T, v {y)jf{w)\V) 



2 6 3 S 



5Af(x)5h^{y)8Af{w) 



(21) 



and the relevant part of the action j7|) that contributes 
to the 3-point function is linear in h pv and quadratic in 
the non-gravitational fields, 



S 



(3) 



d 5 x 
1 



Z -^-{d p <K a - A a p )(d a n a - A%) 



2z 3 



h p<r [-F az F pz + n ^F aa F p p\ 



(22) 



where h p<J is the metric perturbation defined analogously 
to equation ©, viz., g pa = z 2 ^ pa -h pa ). 

To isolate the pion-to-pion elastic stress tensor matrix 
elements from the Fourier transformed 3-point functions 
( x J a (~P2)T pl/ (q)jP(pi)'), we apply the completeness re- 
lation (|15p twice, then multiply by 



P1P2-H 

J U 



(23) 



and take the limit p\ 



and p 2 — *• ml = 0. 



We obtain the transversc-traceless part of the stress 
tensor matrix elements, for T^ifi) at the origin in coor- 
dinate space, 



(7r°( P2 )|f'«'(0)|7r 6 (pi)) 

= 2S ab A^(Q 2 



(24) 



W + Yx («V 



q p q») 



where p = (jp\ +p2)/2 and q = pi — p\. The gravitational 
form factor A v is given by 



(d z ^(z)f v(zf^{zf 



AAQ )= / dzH(Q,z) 



aiftz 



2,3 



ftz 



(25) 



Note that, except for the TC, this form factor is similar to 
the electromagnetic form factor given in [23, [2J] although 
they come from different terms of the action (|7|). H(Q, z) 
is the bulk-to-boundary propagator of the graviton for 
spacelike momentum transfer q 2 = —Q 2 < 0. It is de- 
fined by hnv(q,z) = Ti(Q,z)h^ u {q), where h liV (q, z) is the 
Fourier transform of the metric perturbation hfj, u {x,z) 
In transverse-traceless gauge, q^h^ v - 
linearized Einstein equation becomes 



and h£ = 0, the 



1 



z 3 d z (— d z h fn ,)+q 2 h f , u = Q, 



(26) 



with boundary conditions h(q, e) = 1 and d z h(q, zq) = 0. 
The solution is given by [lj| 



H(Q,z)= l -Q 2 z 2 



KiiQzp) 
h(Qz ) 



h(Qz) + K 2 (Qz)). (27) 



Since TC(0, z) = 1, one can check that ^(0) = 1, which 
is a correct normalization for A„. 

Our procedure obtains the transverse-traceless part of 
the stress tensor; the full stress tensor can have a trace, 
which means there could be a term h(r)^ l/ — q^q v /q 2 )T, 
where T is the trace of T^ v . In general, there are two 
gravitational form factors for spin-0 particles. The ex- 
pression for the pion matrix elements written in terms of 
the two independent form factors is 



(7T a ( P2 )\T^(0)\n b ( Pl )) 



(28) 



c<l li 



24(Q 2 )pV + ^(Q 2 ) (gV" - gV) 



we have calculated A„(Q 2 ), but CV(Q 2 ) = A„(Q 2 )/3 + 
Ctt(Q 2 ) where C is not determined here. 

For oi, the cor resp onding matrix clement is identical 
to the p- meson's [19|. The only difference is that the p- 
meson's wave function tp n is replaced by ipA,m the ai's 
wave function. The A form factor is now given by 



A ai (Q 2 ) = J^H(Q,z)^ 



iipj 



(29) 



The other form factors mirror the p-meson's form factors 
expression. 

Both At,{Q 2 ) and A ai (Q 2 ) are shown in Fig.Q] 



IV. CONSEQUENCES 
A. Radii 



where 



p(z) = 



(d z ^(z)f v(z)H{zf 



2 f2.,2 



%ftz 



2 r 4 



ft* 



(31) 
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FIG. 1: Plot of A with momentum transfer squared in units 
of Aqcd = 1/zo- The red solid line is A v and the blue dashed 
line is A ai 



We find 



(r 2 ) = 0.13(fm) 2 



(0.36 fm) 2 



(32) 



M. Polyakov, in a different model [28|, also found a small 
gravitational radius for the pion, (r 2 ^ = (0.42 fm) 2 . 
The gravitational RMS radius is significantly smaller 
than the electromagnetic radius (r 2 ) c = 0.33 (fm) 2 = 
(0.57fm) 2 , obtained from the AdS/QCD model in chi- 
ral limit J2J]. Compared to the experimental result 
(r 2 ) c = (0.67 fm) 2 |27|, the difference is even more ap- 
parent. This shows, as in the p- meson case, that the en- 
ergy distribution of the pion is concentrated in a smaller 
volume than the charge distribution. 
For a\ , the corresponding radius is 



(r 2 ) 

V ai/gr. 



gra.v 



dz z I 1 



2z 2 



P 



-i,i 



0.15 (fm) 2 = (0.39 fm) 5 



(33) 



As expected, it is smaller than the charge radius 
(r 2 ai ) c = 0.39 (fm) 2 = (0.62 fm) 2 calculated from the 
AdS/QCD model. 



In the limit Qzq <C 1, one can expand 7i{Q 2 ,z) and 
obtain the radius 



(rl). 



-6 



dAn 



dQ 2 

— I dz z" 
4 



1 " ^s ) P(*)> (3°) 



B. High Q 2 limit 

For high Q 2 , the form factor A^ scales as 1/Q 2 . The 



precise expression is given by 

4p(0) 



AAQ 2 ) 



Q 2 



16^ 2 / 2 

Q 2 



(34) 



which follows from the fact that at high Q 2 , the second 
term of the function H{Q, z) in equation l|2?p dominates. 
This term behaves like e~® z . Therefore, one can allow 
Zq — > oo in equation (|25p . and replace p{z) by its value 
at the origin, p(0), and then do the integral analytically. 
One can verify, most easily in the Breit frame, that this 
scaling agrees with the perturbativc QCD prediction. It 
can be shown that [29| 

(pal^T^lpi) ~ Q°, ( P 2\e fiV (q,0)T^\ P i) ~ Q°; 

(35) 
e Aw(<3S 0) gives the helicity-0 component of the spin-2 part 
of the stress tensor, with 

e^M) = -^(2C M («,0)C v (g,0)-C„(?,+)Cv(ff,-) 
-Cm(8, -)&(«,+)). (36) 



Here, £(g, A) is the polarization vector of a spin-1 particle 
of momentum q. Equation (|35p is equivalent to 

Ar.C^l/Q 2 . (37) 

For ai , the high Q 2 behavior of the form factor A ai is 

12|^,„| 2 



a^q 2 ) 



(38) 



Similarly C ai ,D ai ~ 1/Q 4 , while F 0l ~ 1/Q 6 , which 
mirror the scaling results for p-mesons, with the notation 
given in [lj|. 



C. Sum Rules for the GPD 

Deeply virtual Compton scattering process involves a 
target absorbing a virtual photon and subsequently radi- 
ating a real photon. The virtual Compton scattering am- 
plitude can be written in terms of integral involving the 
generalized parton distributions (GPD) H(x,S,,Q 2 ) [30]. 
In a model with quarks, x is the light-cone momentum 
fraction of the struck quark constituent relative to the 
total momentum of the target hadron. 

For spin-0 hadrons, there is only one GPD, defined by 



P +d V ixp+y-/2 I 17/ Vs + , ( Vs\ 

\ Z Z I y+=y ± =0 



■in 



2p + H{x,i,Q 2 ), 



(39) 



where q + = q a + q 3 = — 2£p + . 

There are sum rules connecting this GPD to the grav- 
itational as well as to the electromagnetic form factors. 
The well known sum rule is for the first moment of the 
GPD M 



/: 



where F(q 2 ) is the electromagnetic form factor defined 
by 



(n(P2M q (oh»Mo)\*(pi)) = WHQ 



(41) 



A further sum rule exists because the stress tensor el- 
ement T ++ can be related to the second moment in x 
of the operator whose matrix element defines the GPDs. 
For the pion, the result was given in [201 ] and reads, 



/ dxxH(x,£, Q 2 ) = An 



(Q 2 )-ecAQ 2 )- (42) 



dxH(x,^Q 2 )=F(Q 2 



(40) 



Ref. [2fJ] also uses a chiral Lagrangian to show that 
^4^(0) = 1 (the momentum sum rule) and CV(0) = 1/4. 
One can set £ = in the above equation so that only the 
first term, which is known from the AdS/QCD model, in 
the right-hand side survives. There are also first-moment 
sum rules for the axial vector meson GPDs, whi ch p re- 
cisely parallel the ones given for the p-mesons in [191 ] 



V. CONCLUSIONS 

We have worked out the gravitational form factors of 
pions and of axial- vector mesons using the AdS/CFT cor- 
respondence, and have given the sum rules connecting the 
gravitational form factors, which can also be called stress 
tensor or energy-momentum tensor form factors, to the 
axial sector GPDs. 

A striking numerical result is the smallness of the pion 
radius and of the axial-vector meson radius as obtained 
from A(q 2 ), the gravitational form factor that enters the 
momentum sum rule. This parallels the results for the 
ordinary vector mesons pj| . It suggests that the energy 
that makes up the mass of the meson is well concentrated, 
with the charge measured by the electromagnetic form 
factors spreading more broadly. 

Extensions of the present work include considering fla- 
vor decompositions [j, [y, [32| of the stress tensor and the 
related topic of considering quarks of differing masses, 
to be able to separate contributions with differing inter- 
nal quantum number, including strangeness. Further, we 
would like work on applying the present considerations 
to nucleons JO]. We hope to return to these topics, but 
for now they lie beyond the scope of this paper. 
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